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We show that all pure entangled states of two d- dimensional quantum systems (i.e., two qudits) 
can be generated from an initial separable state via a universal Yang-Baxter matrix if one is assisted 
by local unitary transformations. 
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Entangled state is the cornerstone of quantum informa- 
tion theory [l[ that has many successful applications in 
quantum information processing, such as the revolution- 
ary one-way quantum computer 2] . quantum cryptogra- 
phy 0i dense coding teleportation communica- 
tion protocols and computation Q. Consequently, the 
ability to generate and control quantum entangled states 
has become a far-reaching goal in experimental manipula- 
tion as well as theoretical investigation in recent years. In 
fact, a great number of experiments have been devoted to 
investigating the production of entangled states of pho- 
tons (including the hyperentangled photon pairs) via the 
process of spontaneous parametric down-conversion in a 
nonlinear optical crystal [H , particularly for use in the 
tests of Bell inequalities. 

A fundamental notion in quantum computation (QC) 
is universality: a set of quantum logic gates (i.e., unitary 
matrices) is said to be "universal for QC " if any uni- 
tary matrix can be approximated to arbitrary accuracy 
by a quantum circuit involving only those gates []]]. For 
example, an arbitrary U (2) matrix can be obtained by 
combining the Hadamard gate together with the phase 
shift gates. When such a U(2) matrix is prepared, an ar- 
bitrary state for a single qubit = cos 1 10) +sin |e 1 ^ 1 1) 
can be immediately generated by acting the U (2) matrix 
on the initial state |0). Since entangled states are im- 
portant for quantum information processing, it gives rise 
to a natural question whether an arbitrary pure entan- 
gled state of two qudits (such as the maximally entangled 
state) can be generated via a universal matrix from an 
initial separable state |00)? 

The purpose of this work is to provide a positive answer 
to the above question. We shall study the problem from 
the theoretical point of view based on the Yang-Baxter 
equation, which was originated in solving quantum in- 
tegrable models [j| [lfj , but recently has been shown to 
have a deep connection with topological quantum com- 
putation and entanglement swapping 
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The Yang-Baxter equation, in principle, can be tested in 
terms of quantum optics [1 61 ] - Let us consider a gen- 



eral pure state of two qudits, whick takes of the fol- 
lowing form: |*) = YiijLo (Hj N)a|:?)b, where \i) A and 
\j) b are the orthonormal bases of the Hilbert spaces A 
and B respectively, and /lij 's are complex numbers satis- 
fied the normalization condition X)j i=o I M^i 1 2 — 1- After 
performing an appropriate local unitary transformation, 
the general state | "J) can be recast to a Schmidt-form as 

\ip)2-qudit = ko|00) + m|11) H + K d -i\d -l,d- 1), 

where Kj's (j — 0, 1, d— 1) are the Schmidt coefficients. 
Since \ip)2-qudit is equivalent to |^) up to a local unitary 
transformation, hereafter when refer to generating an ar- 
bitrary pure state of two qudits we mean generating the 
state 1*0)2 qudit- Now we illustrate the problem by start- 
ing from unitary solutions of Yang-Baxter equation in 
the following. 

Unitary solutions of Yang-Baxter equation. The 
Yang-Baxter equation is given by 

R l {x)R l+ i(xy)R l {y) = R i+ i{y)Ri{xy)R i+ i(x). (1) 

Here the notation Ri(x) = Ri t i+i(x) is used, Ri,i+\{x) 
implies Ii ® I2 <8> I3 • • • <8> Ri,i+i (x) (g> • • • <£> I„, Ij represents 
the unit matrix of the j-th particle, and x — e lB is a 
parameter related to the degree of entanglement. Let 
the unitary Yang-Baxter i?-matrix for two qudits be the 
form 



Ri(x) = F(x)[li + G(x) Mi], 



(2) 



where F(x) and G(x) are some functions needed to deter- 
mine later on, 1^ = I, <S)Ij_|_i, and the Hcrmitian matrices 
Mi's (i.e., M{ — At}) satisfy the Hecke algebraic rela- 
tions: (MiM i+1 Mi - M i+ iMiM i+ i)+g(Mi - M i+X ) = 0, 
Mf = aMi + 13 li, with a = d - 2 and /3 = g = 
d — 1. Substituting Eq. (J2]) into Eq. (P), one has 
G(x) + G(y) + aG(x)G(y) = [1 + gG(x)G(y)} G(xy). 
The unitary condition R\{x) = R^ 1 (x) — Riix^ 1 ) yields 
G{x) + G{x- 1 ) + a G{x)G(x- r ) = 0, F(x)F(x^ 1 )[l + 
j3 G{x)G(x^ 1 )] = 0. In addition, the initial condition 
R^x) = It leads to G(x = 1) = 0, F(x = !) = !. As a 
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result, one has 

G(X): 



- F ( x ) 



(d - l)x + x^ 1 



(d- tyx + x- 1 ' ' S ~ J d 
In this work, the d 2 x d 2 matrix M is realized as 

d-1 d-1 d-1 

M = ^P r ®P r = \*3)^,7+^l (3) 



r=l 

where i + r = Mod \i 



iJ—O r—1 

■ r, d] , and 



d-1 

£ 

i=0 



0,1,. 



(4) 



are the circulation matrices that transform the basis 
{|r), |r+l), |r + 2), |2), |1), |0)} of a qudit to the basis 
{|0), |1), |2), \d- 3), |d- 2), |d - 1)}. The operator P r 
can be realized through the multiplication of permuta- 
tion operators P fc , fe+ i = (I - \k)(k\ - \k + l)(fc + 1|) + 
\k)(k + 1\ + \k + l)(k\ of a single qudit, for example, 
Pi = P(2-2,d-i"-IP2,3lPi,2- Moreover, the traceless matri- 
ces P r 's satisfy the following interesting relations: 



P-r — (-Pi) 7 PnPi — PiPn — P]VIo 



d[m+n,d] ■ 



Let Pq = Y2i=o denote the d x d unit matrix, we 
eventually arrive at the unitary Yang-Baxter matrix for 
two qudits as 

Ri{x) = -{[(<*- l^ + aT 1 ]P <g>P 



d-1 

-(x-x- 1 ) VP,. 



r=l 



P, 



(5) 



which has not been reported in the literature. Our main 
result of connecting unitary Yang-Baxter matrices with 
entangled states of two qudits is the following Theorem. 

Theorem: All pure two-qudit entangled states 
\'4 > )2-qudu can be generated from an initial separable 
state 1 00) via a universal Yang-Baxter matrix R(x) if 
one is assisted by local unitary transformations Ua®Ub 
and Va®Vb, namely, 

= [V A ® V B ] R(x) [U a ® U B ] |00). (6) 

Here the local unitary transformation Va ® Vb is intro- 
duced in order to transform a two-qudit state into its 
Schmidt-form. 

Proof. We would like to provide analytical proof for 
the case with d = 2 and numerical proof for the cases 
with d = 3 and 4. 

i) For d = 2, in this case Pi = ( ^ ^ . When one 

acts R(x) directly on the separable state |00), he then 
generates the following family of states 



)yb = - 



(x + a; -1 )|00) - (x-x-^lll) 



(7) 



In Ref. 171 ] . the generalized concurrence (or the degree 
of entanglement [181 ]) for two qudits is given by 



C = 



l-/i 



(8) 



■ + \Kd- 



where/i = Tr[p^] = Tr[p|] = M 4 + |«i| 4 
Pa and ps are the reduced density matrices for the sub- 
systems. For d = 2, one easily has C = 2|k ki|. Obvi- 
ously, \4>)yb has already been in the form of \ip)2-qubit = 
Ko|00) + ki|11), with k = (x + x~ x )j2 = cos# and K\ = 
~(x — x~ 1 )/2 = —i sin 8. The degree of entanglement for 
the state \i/))yb equals to C = 2\k k 1 \ = \ sin(29)\, which 
may range from to 1. Thus, for the case of two qubits, 
all pure states can be generated from 1 00) directly via a 
universal Yang-Baxter matrix R(x). By the way, when 
8 = 7r/4, the state \iP)yb becomes the maximally entan- 
gled state, or the Bell state \ip) Be ii = (l/v / 2)(|00>— i]ll>). 
ii) For d = 3. in this case 




When matrix R(x) is acted directly on |00), it yields the 
following family of states 

\1>) YB = \[{2x + x-^OO) -(x- x-^dll) + |22))], (9) 
whose generalized concurrence reads 




1 \2x- 
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When |2a; + .T _1 | = \x— x~ x \, namely x — e lTr / 3 , the state 
\iP)yb becomes the maximally entangled state (here we 
would like to call it as the GHZ state) of two qutrits as 
Wighz = 7§M°°) + I 11 ) + I 22 )]- In general, if one acts 
the unitary Yang-Baxter matrix R(x = e l7T ^ 3 ) on the ba- 
sis {|00), |01), |02), 1 10), |11), |12), |20), |21), |22)}, he will 
generate nine complete and orthogonal maximally entan- 
gled states of two qutrits. 

It is easy to check that the generalized concurrence C 
ranges from to 1 when the parameter 8 runs from 
to it. However, this fact does not mean that \iP)yb is 
an arbitrary state of two qutrits, because \ip)2-qutrit has 
at least two free parameters while \iP)yb contains only 
one. Actually, the entanglement property of a two-qutrit 
system is completely charecterized by two entanglement 
invariants h = Tr[p 2 A ] = Tr[p|] = |k | 4 + |«i| 4 + M 4 
and h = Tr[p 3 A ] = Tr[p B ] = \k \ 6 
equivalently, 



«i 



K2 



or 



Ii 



o -g(l-/ 2 ), 



(10) 



FIG. 1: (Color online) In the /( — -Z2 coordinate, the separable 
states, such as |00), locates at the origin O = (0, 0); the max- 
imally entangled state locates at the point G = (1, 1), which 
is the point farthest from the origin; and the entangled state, 
such as 775(100) + [11)), locates at the point B = (|, §|). 
One may plot the points (7i,/ 2 ) for the states \iP)'yb = 
R(x) [\0}a ® (cos^|0)s +B3Htp\l) B )] = R(x) [U A <S>U B ] |00), 
and they perfectly recover all the red region of figure. 



FIG. 2: (Color online) In the I[ — I' 2 — l' 3 coordinate, we have 
plots points 12,13) of the two-qudit state \ip)2-qudit for 
d = 4, see the blue region. By randomly taking 10 7 val- 
ues of 8, ifi and ip2, one may also plot points I^I'z) 
for the states \ip)' YB = R(x) \$}a ® |$)s = R(x) |0)a ® 
(cos<^i|0)b + sin ipi cos ^2 1 l)s + +sin^i sini/52|2)s), which 
perfectly recover all the blue region of figure. 



where the normalized entanglement invariants I[ , I 2 £ 
[0,1]. 

In Fig.l, we have plots points (I[,I 2 ) f° r the two- 
qutrit state \il))i- qu trit = ko|00) + Kijll) + k 2 |22) by 
randomly taking 10 7 values of kq, k±, and K2, see the 
red region of figure, whose contour lines form a curved 
triangle AOBG. One may observe that for a fixed value 
of I[, there are different values for I 2 , therefore I[ is 
not enough for characterizing the entanglement property 
of two qutrits. In the /{ — I 2 coordinate, the separa- 
ble states, such as |00), locate at the origin O = (0,0). 
The maximally entangled states (or say the GHZ states), 
such as Wighz = TjdOO) + |11) + |22)), locate at the 
point G — (1,1). And the entangled states, such as 
-7=(|00) + |11)), locate at the point B = (f,§). The 

contour line OB corresponds to the states \ip) 2 _ qutr it = 
cos£|00) -I- sin £| 11), the point I£) runs from O to B 
when £ runs from to 7r/2; The contour lines OG and 
GB correspond to the states \ip} 2 - qu trit = cos£|00) + 
sin f (| 11) + |22>)/ v / 2, the point (I{J^) runs from O to 
G when £ runs from to 7r/3, and from G to B when £ 
runs from 7r/3 to 7r/2. 

The state \ip)yB is one part of the state IVOl-gutrtf 
when 9 runs from to it/2, the point (J{, 1%) runs from 
O to G, then runs along the line GB towards to point B 
and finally stops at a point (|, ||), which corresponds to 
the state \ip) Y B = |(|00) - 2|11) - 2|22>). Namely, when 
one acts R(x) directly on the state |0)^|0)s, he cannot 
get all pure state of two qutrits. However, numerical 
computation shows that if one acts R(x) on the state 
\0)a <8> (cos(p\0)b + sin (p\1)b), he can indeed obtain all 



pure two-qutrit states: by randomly taking 10 7 values of 
9 and ip, one may plot points (/{, I 2 ) for the states 

W)'yb = &( x ) [|0)a® (cos ip\0) B + sin <p\l) B )], 

which perfectly recover all the red region of Fig. 1. This 
means that if one is assisted by the local unitary trans- 
formations Ua ® Ub and Va <8> Vg, with Ua = I and 
(7 B =cos^(|0)(0|-|l)(l|) + sin^(|0)(l| + |l)(0|) + |2)(2|, 
he then can generated all pure two-qutrit entangled states 
in the following way: \ip)2-qutrit = [Va®V b ] R[x) [U a ® 
Ub] 1 00). This ends the numerical proof for the case with 
<2=3. 

iii) For d > 4, there are d — 1 normalized entangle- 
ment invariants for the two-qudit state \ip) 2—qutrit, i.e., 
I'i = J^ii 1 ~ hi with Ij = Tr[^ +1 ] = Tr^ 1 ] = 
£r=oK| 2(j ' +1) , (j = 1,2,..., d- 1). When acting on 
the separable state |00), the unitary Yang-Baxter ma- 
trix R(x) generates the following family of states 

d-l 

W)YB = ~ d {[{d - 1)X + -(X- X' 1 ) ]T |ji)}. 

J = l 

If one requires the state \ip)yB to be the maximally en- 
tangled state (or the GHZ state), he must set \(d— l)x + 
x~ 1 \ = \x — £ |, namely, 

cos(20) =1-5, 

For d — 2, d — 3, and d = 4, one has 9 = ir/4, ir/3, 
and 7r/2, respectively. However, the above condition is 
not valid for d > 5, because one will have |cos(20)| > 



4 



1 when d > 5. This fact implies that the maximally 
entangled two-qudit states can be generated when we act 
R(x) directly on the separable state 1 00) for d < 4. 

Similarly, numerical results show that all pure two- 
qudit entangled states can be generated in the following 
way: \ip) 2 - q udit = [V A <8> V B ] R(x) [U A <8> U B ] |00) = 
[V A ®V B ] R(x) \*)a®\*)b, where \$} A = U A \0) A = \0) A 
and = U B \0) B = cos^ilO)^ + sini^i cos^ 2 |1)b + 

• ■ • + sinyi sin(/j2 ■ • • simpd—2\d — 2)b- In particular, for 
d — 4, one has |$) s = cos^ilO) b + sin^i costal) s + 
+ sintpi sin^>2|2) s- Numerical proof of the Theorem for 
d = 4 is provided in Fig. 2. In Fig. 2, we have plots 
points I2, 13) for the two-qudit state IV^a-gadit = 
Ko|00) + «i|ll) + K2|22) + +K3|33) by randomly taking 
10 7 values of Kj's, see the blue region of figure. By ran- 
domly taking 10 7 values of 6, ipi and tp 2 , one may also 
plot points (I[,l2,I^) for the states 



Yyb 



R(x) \§)a®\$)b, 



(11) 



which perfectly recover all the blue region of Fig. 2. 

In conclusion, we have shown that all pure entangled 
states of two qudits can be generated from an initial sep- 
arable state 1 00) via a universal Yang-Baxter matrix if 
one is assisted by local unitary transformations. Even- 
tually, we would like to point out that the spirit of a 
unitary matrix assisted by local unitary transformations 
as shown in Eq. ([6]) or Eq. (| 1 1 1) coincides with the spirit 
of entangling power, which is a quantitative measure how 
much entanglement capability a given unitary operator 
has in the context of quantum information. The concept 
of entangling power is first introduced in Refs. OHil, 
which is defined as 



e v {U) :=E(U \$} A ® |*) B ) 



(12) 



where E(\^f)) = 1 — Trfp 2 ^], the overbar stands for the 
average over all the product states, and it can be sim- 
plified as e p (U) = (d/d + l) 2 [E(U) + E{US) - E{S)}, 
with S = JY j^_ Q is the permutation operator 

of two qudits. The entangling power has been useful 
for the stud y o f quantum evolutions and Hamiltonians 
[H 0, ED, lH III Hi|, and been also applied to 
some quantum chaotic systems [2(| 27, [2^, [29}. Actu- 
ally, E(\^)) = 1 — Trlp 2 ^] is the entanglement invari- 
ant J{ up to a normalized constant d/(d — 1). Simi- 
larly, based on the entanglement invariants of two qu- 
dits one may define a series of entangling powers as 
e^(U) := E 3 {U with E, 

will investigate subsequently. 

This work is supported in part by NSF of China 
(Grants No. 10575053 and No. 10605013), Program for 
New Century Excellent Talents in University, and The 
Project-sponsored by SRF for ROCS, SEM. 



I'j, which we 



[1 



[6 



* Electronic address: chenjl@nankai.edu.cn 
t Electronic address: gemlQ nankai . ed u . cn | 
M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2000). 

R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 
C. H. Bennett, and S. J. Wiesner, Phys. Rev. Lett. 69, 
2881 (1992). 

C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. 
Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 
(1993). 

R. Prevedel, P. Walther, F. Tiefenbacher, P. Bohi, R. 
Kaltenbaek, T. Jennewein, and A. Zeilinger, Nature 445, 
65 (2007). 

[7] P. G. Kwiat et al, Phys. Rev. Lett. 75, 4337 (1995). 
[8] J. T. Barreiro, N. K. Langford, N. A. Peters, and P. G. 

Kwiat, Phys. Rev. Lett. 95, 260501 (2005) 
[9] C.N. Yang, Phys. Rev. Lett. 19, 1312 (1967). 
[10] R. J. Baxter, Exactly Solved Models in Statistical Me- 
chanics (Academic Press, London, 1982). 
[11] A. Y. Kitaev, Annals Phys. 303, 2 (2003) 
[12] M. H. Freedman, A. Kitaev, and Z. Wang, Commun. 
Math. Phys. 227, 587 (2002); M. H. Freedman, M. 
Larsen, and Z. Wang, Commun. Math. Phys. 227, 605 
(2002). 

[13] J. Preskill, Lecture Notes on 

Topological Quantum Computation 

( http: //www. theory, caltech. edu~preskill/ph219 ). 

[14] L. H. Kaufhnan and S. J. Lomonaco Jr., New J. Phys. 6, 
134 (2004). 

[15] J. L. Chen, K. Xue, and M. L. Ge, Phys. Rev. A 76, 
042324 (2007); J. L. Chen, K. Xue, and M. L. Ge, Ann. 
Phys. 323, 2614 (2008). 

[16] S. W Hu, K. Xue, and M. L Ge, Phys. Rev. A 78, 022319 
(2008). 

[17] S. Albeverio and S. M. Fei, J. Opt. B: Quantum Semi- 
class. Opt. 3, 223 (2001). 
[18] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 

(1997) ; W. K. Wootters, Phys. Rev. Lett. 80, 2245 

(1998) . 

[19] P. Zanardi, C. Zalka, and L. Faoro, Phys. Rev. A 62, 

030301 (R) (2000). 
[20] P. Zanardi, Phys. Rev. A 63, 040304(R) (2001). 
[21] W. Dur, G. Vidal, J. I. Cirac, N. Linden, and S. Popescu, 

Phys. Rev. Lett. 87, 137901 (2001); J. I. Cirac, W. Dur, 

B. Kraus, and M. Lewenstein, Phys. Rev. Lett. 86, 544 

(2001) ; G. Vidal, K. Hammerer, and J. I. Cirac, Phys. 
Rev. Lett. 88, 237902 (2002); W. Din, G. Vidal, and J. 
I. Cirac, Phys. Rev. Lett. 89, 057901 (2002). 

[22] M. S. Leifer, L. Henderson, and N. Linden, Phys. Rev. A 

67, 012306 (2003). 
[23] X. Wang and P. Zanardi, Phys. Rev. A, 66, 044303 

(2002) . 

[24] A. M. Childs, D. W. Leung, F. Verstraete, and G. Vidal, 

Quantum Inf. Comput. 3, 97 (2003). 
[25] B. Kraus and J. I. Cirac, Phys. Rev. A 63, 062309 (2001). 
[26] R. Demkowicz-Dobrzahski and M. Kus, Phys. Rev. E 70, 

066216 (2004). 

[27] A. J. Scott and C. M. Caves, J. Phys. A 36, 9553 (2003). 



■5 



[28] R. O. Vallejos, P. R. del Santoro, and A. M. O. de (2006). 

Almeida, J. Phys. A: Math. Gen. 39, 5163 (2006). 
[29] R. F. Abreu and R. O. Vallejos, Phys. Rev. A 73, 052327 



